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t INTRODUCTION 

The  years  following  the  Second  World  War  have  seen  a substantial  development  in 
understanding  the  search  and  detection  game  involving  two  military  units  searching  for 
each  other  with  some  specified  detection  hardware.  It  has  become  apparent  that  Koopman’s 
attempt  - a good  part  of  which  stands  canonized  in  his  1946  book  - establishes  a suitable 
framework  for  the  description  of  these  operations.  Much  of  the  early  work,  however,  fails 
to  adequately  accommodate  the  modified  rules  of  the  game  that  technological  development  as 
well  as  increased  analytical  exigencies  have  since  imposed  upon  us. 

It  is  the  contention  of  the  author  that  notwithstanding  the  proliferation  of  papers  that 
have  been  written  on  the  subject  attempting  to  account  for  the  complexities  of  modern  search 
and  detection,  a clear  framework  malleable  enough  to  handle  a large  variety  of  outstanding 
real  life  problems  has  failed  to  emerge.  Therefore,  the  purpose  of  this  paper  will  be  to 
revisit,  criticize,  slightly  modify,  and  subsequently  expand  the  original  Koopman  framework . 

Perhaps  it  is  worth  noticing  from  the  onset  that  we  have  chosen  to  address  little,  if 
anything,  of  this  work  to  the  question  of  comparison  with  the  state  of  the  art.  W'e  have  done 
so  because  we  are  more  interested  in  presenting  an  independent  framework,  who^e  right  for 
survival  we  shall  heretofore  attempt  to  argue,  rather  than  improve  on  limited  treatments 
available  in  the  literature. 

1.  THE  KOOPMAN  APPROACH  REVISITED 

In  the  original  Koopman  formulation,  the  game  of  sea-ch  and  detection  includes  two 
players,  the  searcher  and  a target.  The  two  players  engage  in  some  sort  of  geometrical 
search  pattern,  whereby  each  one  of  them  moves  along  deterministic  trajectories.  This 
kinematical  part  of  the  game  is  well  described  mathematically  by  the  relative  track  of  the 
target, 


K = S(t)  ; n = n(t) 


(l-D 


While  moving  on  its  trajectory  the  searcher  is  assumed  to  be  looking,  either  through 
glimpses  or  continuously,  for  the  target.  The  work  'Tooking"  is  used  here  to  mean  attempt 
of  detection  via  some  hardware,  be  it  visual,  radar,  sonar,  or  any  other  means  whatso- 
ever. This  part  of  the  game  is  pure  / probabilistic  in  nature  and  the  central  mathematical 
entity  to  be  considered  is  the  instantaneous  probability  density  of  detection  y(t)  , where: 


Y ( t) 


Y 


(t) 


+ n2(t) 


(1-2) 
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Notice  that  in  equation  ( I -2),  v{t ) is  not  taken  to  depend  on  the  relative  bearing  of  the 
target.  While  not  a necessary  assumption,  it  certainly  is  a very  reasonable  one  for  most 
deter  ion  hardware. 


Once  v(t)  is  given,  the  probability  of  detection  for  continuous  looking  through  time 
t is  given  by: 


- [ dxy  (t) 


p(t) 


d-3) 


In  principle,  equation  (1-3)  provides  the  probability  of  detection  for  relative  tracks  of  any 
given  complexity.  In  practice,  one  is  limited,  of  course,  by  one's  capability  of  estimating 
the  exponent  in  (1-3)  along  the  track.  In  fact,  the  task  of  evaluating  the  integral  in  question 
is  more  often  than  not  a quite  formidable  one,  and  correspondingly,  solutions  are  difficult 
to  get  for  all  but  a few  instances  which,  in  virtue  of  their  simplicity,  are  hardly  realistic. 
It  is  also  evident,  that  a solution  to  equation  (1-3)  is  less  than  realistic  in  that  it  fails  to 
take  into  account  the  random  nafre  of  target  arrival. 


A somewhat  more  significant  result  obtains  if  one  is  willing  to  restrict  the  Koopman 
treatment  by  making  the  following  strong  assumptions: 

• The  target's  position  is  uniformly  distributed  over  some  region  of  size  A in 
the  ocean. 

• The  searcher's  path  is  random  in  A in  the  sense  that  it  can  be  thought  of  as 
having  its  different  (not  too  near)  portions  placed  independently  of  oue  another 
in  A. 


Then,  the  probability  of  detection  is  given  by  the  Koopman  random  search  formula, 

- WL/A 

p(I.)  =l-o  , (1-4) 


where  \Y  stands  for  the  effective  search  width  of  the  detection  law  and  L represents 
the  total  length  of  the  searcher's  path  within  A . Let  us  emphasize,  however,  that  while 
this  development  does  indeed  eliminate  the  difficulty  of  integrating  v(t)  along  the  relative 
track  of  the  target,  it  only  does  so  at  the  expense  cf  sacrificing  che  realistic  and  eften 
crucial  concept  of  correlation  along  the  path  of  the  searcher. 

2.  THE  KOOPMAN  APPROACH  MODIFIED 

Tiie  Koopman  framework  lias  proved  quite  successful  in  a large  variety  of  naval  war- 
fare problems,  and  has  served  as  a guideline  to  many,  more  sophisticated  later  treatments. 
It  is  within  the  nature  of  our  stated  purpose,  however,  to  investigate  the  shortcomings  of 


the  formulation  rather  than  emphasize  its  virtues.  Coirespondingly,  we  shall  concentrate 
in  the  following  upon  tl.e  major  limitations  of  the  Koopman  framework. 

As  will  be  recalled,  the  Koopman  theory  takes  the  a posteriori  view  of  the  search  and 
detection  problem  whereby  both  target  and  searcher  trajectories  are  assumed  known.  More 
often  than  not,  however,  such  is  nrt  the  case , In  fact,  the  path  of  a military  unit  is  well 
’■epresented  by  a juxtaposition  of  deterministic  pieces  of  path,  the  process  of  passing  from 
one  tc  any  other  in  time  ueing  a r adorn  one  to  be  sampled  from  some  given  distribution 
dependent  on  both  tactical  and  physical  parameters.  The  question  to  be  answered,  then, 
is  one  concerning  the  manner  in  which  the  probabilistic  description  inherent  to  the  arrival 
process  is  to  be  incorporated  in  the  Koopman  equation.  While  we  do  recognize  that  an 
an  attempt  at  addressing  this  question  exists  in  the  Koopman  theory,  it  is  important  to  note 
that  limitations  in  the  mathematical  formulation  thereof  do  not  allow  the  preservation  of 
correlations  along  the  path  of  the  searcher. 

It  is  also  worth  noticing  thit  the  two  parts  of  the  game  as  envisioned  in  the  Koopman 
formalism  are  not  really  independent  of  each  other,  in  that  looking  can  produce  feedback 
on  the  trajectories  of  both  players.  In  this  sense,  the  framework  is  void  of  dynamics. 

The  formulation  that  we  propose  and  develop  to  some  extent  in  this  work  is  designed 
to  hopefully  handle  each  and  all  of  the  points  raised  against  the  Koopman  framework.  It 
essentially  consists  in  determining  various  physical  quantities  relevant  to  the  search  and 
detection  game  f<  * any  of  a family  of  simple  target  paths,  each  one  of  which  is  assigned 
a given  probability  of  havi.g  been  chosen  from  the  family.  The  actual  value  of  the  physical 
quantity  of  interest  is  then  obtained  by  averaging  over  the  family  parameters. 

The  main  conceptuaT  virtue  ct  this  framework  <s  that  it  is  more  responsive  to  the 
actual  restrictions  that  realit*  imposes  upon  the  search  and  detection  game.  This  added 
responsiveness  gets  reflected  not  only  in  the  possibility  of  choosing  the  path  probabilities 
to  suit  the  real  case  at  hand,  but  also  in  the  freedom  A choosing  the  family  in  such  a man- 
ner that  it  be  indicative  of  the  operational  features  of  the  actual  setup  while  still  simple 
enough  to  allow,  in  principal  at  least,  for  tractable  solutions  to  equation  (1-3).  It  is  inter- 
esting to  note  that,  unlike  the  Koopman  formulation,  ail  of  these  conceptual  advr  ges  are 
available  without  the  loss  of  correlation  along  the  path  of  the  searcher.  In  fact,  as  we 
shall  discuss  later  on  in  some  detail,  most  of  the  relevant  correlations  are  properly  ".aken 
into  account  by  the  theory.  Statements  co.icerning  the  less  conceptual  advartages  the  frame- 
work offer , are  relegated  to  conclusions  when  the  reader  shall  be  more  familiar  with 
the  lanwiage. 

3.  OUTLINE 

For  reasons  of  simplicity  we  have  chosen  to  break  the  presentation  to  follow  up  into 
two  parts,  in  the  first  chapter  we  shall  address  the  relatively  independent  question  of 
determining  realistic  probabilities  of  arrival,  both  because  they  are  relevant  in  them- 
selves and  because,  as  it  turns  out,  they  are  the  necessary  building  blocks  for  calculating 

-3- 


probabilities  of  detection.  The  second  chapter  of  this  paper  will  shew  how  nonuniform  di.; 
tributions  of  target  arrival  cau  be  incorporated  into  the  detection  problem.  Since  we  suc  - 
ceed in  doing  so  without  breaking  up  the  searcher  path  in  independent  pieces,  correlations 
along  the  path  of  the  searcher  are  fully  preserved  in  the  formalism. 


-4- 


— w u£$ 


PART  I.  NONUNIFORM  PROBABILITIES  OF  ARRIVAL 


A.  THE  GENERAL  FRAMEWORK 

We  pursue  the  progiam  outlined  in  the  Introduction  to  this  pa  pc-  by  addressing  Lrst 
the  relatively  independent  question  of  calculating  realistic  probability  densities  of  target 
arrival.  We  thus  let, 

F(r,u>)  = 0 (A-l) 

be  an  arbitrary  family  of  curves  parameterized  by  the  set  to  and  define  the  probability 
a(to)  that  the  target  has  chosen  to  approach  the  searching  area  along  some  path  to  . 
Naturally, 


/ dw  a (to)  = 1 


(A -2) 


Now,  if  (p  , 1 1 pQ  , t ) represents  t he  probability  density  of  arriving  at  ( p , t ) 
from  (p  , t^)  along  path  <o  , we  write  for  the  probability  density  of  arrival  at  (p  , t) 
having  started  at  (p^  , t ^ ) , 


CCp.t|p0,t0)  - j du)  a ( to ) G ( p , 1 1 p ^ , t q ) 
To  proceed,  we  now  specify 

Ga)  C P » T I » 1 o ) = fipU,t)-p) 


(A -3) 


(A -4) 


where  r f to , t)  describes  the  position  of  the  moving  target  along  the  path  to  . Corres- 
pondingly, 


In  fact, 


G(P,t  = j dto  d[wj  j . 


J dto  a(ti)  r(u>,t)  = / dto  a(u>)  j dp  p 6^rfw,t)-p  j 


(A -5) 


- /dp  p J dto  a (to)  6 (to,  t) -p  j , 


-5 


which  identifies  G (p  . 1 1 pQ  , t ) given  inequation  (A-5).  For  the  reader  who  is  not  im- 
pressed with  the  rigor  of  our  derivation,  we  note  that  one  can  easily  continue  the  argument 
to  show  that  (A-5)  also  has  the  appropriate  higher  moment. 

Using  the  properties  of  the  Dirac  delta  function,  on  ■ immediately  shows  that 
G(p,t  | is  properly  normalized  at  all  times  if  a(to)  satisfies  equation  (A-2).  Indeed, 

using  equation  (A-5),  we  have: 

fd$  G , 1 1 pQ  » t Q ) = J doj  a(<a) 

and  hence, 

fdp  G(p,t|p  ,tQ)  = 1 . 


It  is  noteworthy  that  the  global  correlation  expressed  in  the  requirement  that  the  target 
reach  (p  , t)  having  started  at  ( ) is  here  automatically  translated  into  purely  geo- 

metrical restrictions  to  be  imposed  upon  the  family, 

?(w,t0)  = pc 

as  well  as  upon  the  integration  in  equation  (A-5).  This  interesting  feature  of  the  formalism 
shall  prove  of  some  importance  later  on  in  that  it  offers  a convenient  way  of  including  cor- 
relations into  the  framework. 


B.  THE  CASE  JF  THE  ISOTOPIC  WALK 

To  start  the  framework  off,  a sensiin  , choice  G,-  the  family  of  paths  is  made.  Fcr 
purposes  of  illustration  and  because  they  are  indie.  _xvc  of  many  re  listic  cases,  we  analyze 
the  family  of  all  r:  lorn  walks  starting  at  the  origin  and  parameterized  by 


it  \ i 

w=|n,(0,;v  ;i.J  j . m . . 

[ V1  1 Vi-1.2 n!  (B-‘> 

Here,  as  displayed  in  figure  1,  the  set  (0  . ,v  - ,a  . ) stands  for  the  angle,  speed  along,  and 


time  interval  of  the  ith  step  m the  walk,  while  n represents  the  number  of  steps  con- 
sidered. Specifically,  for  a random  walk  of  n steps  the  family  equation  reads: 
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t 


f 


n 

I 

i = l 


a-  cos  8. 
i 1 


n 

y (t)  = I a.  sin  6 
n i = l 3 1 


(B-2) 


a . 

i 


vi' 


t I t.  . (B-3) 

i = l 


In  this  section  we  specialize  t the  simple  case 

= a ; = v ; i = l,2,...,n 

and  correspondingly  measure  time  in  units  of  t = ct/v  , 

t = nx 


(B-4) 


(B-5). 


l ie;,  i:  Tin;  family  or  ranidm  ivai.ks 


Therefore,  the  associated  quantity  a(oo)  becomes  the  joint  probability  density 

,...,6  ) that  the  family  parameters  take  on  the  set  of  values  j 

i 2 n 1 2 n 

respectively,  and 
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(B-6) 


2tt  2tt 

“o  1/ 

0 0 


/ 11 

6(1  a cos  0 . 

\i  = l 

for  (q,  0)  the  cartesian  components  of  P 

We  now  make  the  simplifying  assumption  that  the  decision  of  choosing  any  particular 

value  for  some  given  0 is  independent  of  the  past, 

i 

a(01,02,  • ■ • ,en)  = a1(e1)a2(02)  ...  an(Bn)  . (B-7) 

Because  (B  -7)  is  tantamount  to  assuming  that  the  probability  distribution  for  any  given 
angle  6 . is  independent  of  the  point  in  the  ocean  at  which  the  walking  unit  finds  itself,  we 

shall  refer  to  (B-7)  as  the  "isotopic  walk  case."  We  hasten  to  add,  however,  that  while 
of  some  methodological  importance  here,  the  isotopic  walk  assumption  is  by  no  means 
necessary.  In  fact,  realistic  correlations  along  the  path  of  the  walk  will  be  properly  taken 
into  account  in  a later  chapter  of  this  paper. 

Finally,  to  make  matters  as  simple  as  possible,  we  specify; 

ai(ei)  = 1/2  -IT  ; i = l,2,.  . . ,n  (B-8) 

while  mentioning  that  the  minor  restriction  imposed  by  (B-8)  can  be  lifted  with  no  difficulty 
at  all. 

Hence, 


G(p, t| 0,0)  = J 


G('p,t|0,0) 


/n 

61  Z a cos 


.*  11 
u 


a sin  0.- 
1 


(B-9) 
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f 
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To  cast  equation  (B-9)  into  a more  manageable  form,  we  make  use  of  the  spectral  repre- 
sentation for  the  Dirac  delta  function, 


6 (x) 


2n 


J dk  elkX 


and  have: 


(B-10) 


(\n+?  “ oo 

jA  ' f d Ce'1?0  J dne'1^ 


(B-ll) 


/ 


2tt 


dGe 


i£ci  cos  e in«  s in  0 


n 


Recognizing  that 
2tr 


f , i^o  cos  0 sin  0 

_/  d0e  e 


= 2tt  dQ^a 


/" 2 
/£  + 


C + n 


(B-12) 


where 


J0(o/d  * 1!) 

G(p,t|?,0)  = (—)  f dCe'l£,°  f dne 
\2ttJ  j 

\ / -00  -oo 


denotes  the  zeroth  Bessel  function  of  the  corresponding  argument, 

nrt  rr» 

2 


JQn  l ol/k  + r n 


(B— 13) 


We  shall  find  it  convenient  to  pass  to  polar  coordinates  in  ( E, , n ) , 

5 - x cos  <J>  ; n = x sin  4> 

and  explicitly  perform  the  angular  integration, 

G(p;t|(JfO)  = ~ J dx  xJo^x,42  + ^ Jon(ax)  . 


(B-14) 
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An  immediate  off  shot  of  (B-14)  is  that 


G(p,t|0,0)  - 0 if  /o  + p~  > na  (B-15) 

as  expected. 

Wc  can  now  start  the  simple  but  lengthy  job  of  calculating  G ( p , 1 1 0 , 0 ) for  various 
values  of  n , Thus,  for  n=2,  (see  figure  (2), 


PIG.  2:  GliOMHTRY  GONCPRNING  TUI:  ARRIVAL 

PROBABILITY  FOR  »=? 


G(p, 2x1(5, 0)  = — — 
4 xr 2 A 


(B -16) 


where  A represents  the  area  under  the  triangle  formed  by  the  triplet  + P ;a;aj 

Naturally,  in  accordance  with  equation  (B-15),  the  probability  density  of  arrival  for  n=2 
vanishes  whenever  the  triplet  fails  to  form  a triangle.  We  also  note  that  for  the  marginal 


case  of  measure  zero. 


A 


2 


2 

+ p 


2a 


where  A-0  , the  probability  densit"  of  arrival  becomes  infinite. 
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To  handle  the  case  n >2  , we  could  make  repeated  use  of  the  equality 


f'*»  Jo( 
0 ' 


Jz~ 


£ - lit,  cos  6 j = JqCZ)Jq(?) 


(B-17) 


to  reduce  the  power  of  the  Bessel  function  to  n-2  , and  then  employ  (B-16),  but  restrain 
from  doing  so  here.  Rather,  we  develop  in  the  next  section  a more  general  technique  of 
handling  the  calculational  problem  that  shall  provide  answers  at  a much  smaller  expense, 
and  at  the  same  time,  indicate  the  manner  in  which  generalizations  to  the  simple  example 
considered  here  are  to  be  approached. 

C.  THE  EQUIVALENT  ITERATIVE  FORMULATION 

To  begin  with,  recall  the  result  (B-13)  from  the  previous  section, 


>(p»nt  j 0 , 0)  = f d£e  dne‘irip  JQn 


(C-l) 


but  avoid  explicit  integration  on  ( £ , n)  . Instead,  make  use  of  the  convolution  theorem 
for  the  Fourier  transform, 


f 


dxe 


ikx 


f(x)g(xl  = ± J dZ  f(£)g(k-*) 


where 


-=  / 


dxe  f(x) 


and 


g(£) 


dxe  ^ ^ X g(x) 


onto  the  product 


f /2  2 

a/f;  +p 


to  have: 


•11 


, _►  - 00  . 00  , 2 2 

G(p>nT | 0 , 0)  = — L__  f dfe'1?0  f Hv  6 l<*  -Cp-y)  ] 

(2tt)3  J J n- 7 


'a  -(p-y) 


(t/J 


' (P_y)  ) f dve  iVy  J, 


'5  +v 


and  continuing, 


G(p,m|0,0)  - -J—  7dv  (p-y);]  r ^ - ivy  , 

<2">3i  .n—i'L 


(2tr)3 


2/a  -(p-y) 


J dx  Sfa-x-A  (p-y)2  j+  s/a-xWa-ip-y) 


j dye"llJX  J^'1/  a^V 


Hence,  recognizing  that 


,2  2 2 
<5  (a-x)  + (p-y)  -a 


2 /a  "(p-y) 


— 6(a-x- /o.7  - (p-y)  2 


+ 6 (a-x  + /a  -(p-y) 


2 2 

and  correspondingly  dropping  the  step  function  0[a  -(p-y)  ] whose  value  is  always  one 
on  the  circle 


2 2 2 
(o-x)  + (P-y)  = a 


the  probability  density  of  arrival  at  p in  n steps  becomes: 


G(p,nr|0,0)  - f dxdy  6[(o-x)  + (o-y)  -a  | * 

f2 L I J 

] dpe-i'131  j dve'1^  j . 


Finally,  since  in  virtue  of  equation  (C-l), 


j due' J dve~iVy  J 


n -l(  / 

) l a/lJ 


2 2 
+ V 


= C 2tt)  G[t,  (n-1)  T I?,  0 ) , 


for  r the  two-dimensional  vector  of  components  ( x,  y)  , our  main  result  reads: 


G(p»nx |M) 


(C-3) 


A quick  glance  at  (C-3)  will  now  convince  the  reader  that  the  physical  content  em- 
bodied therein  is  surprisingly  simple:  the  probability  G(p ,dt  | ft,  0)jdp  of  arriving  at  ~p 
in  n steps  having  started  at  the  origin  is  given  by  the  probability  G (r  , (n- 1) r | ft, 0 j dr  of 
arriving  at  r in  (n-1)  steps  having  started  at  the  origin  folded  into  the  probability 


.1  •+■  -*•  2 2 -y  — _y 

--  6[(p-r)  -a  ] dp  of  taking  the  last  step  from  r to  p and  integrated  over  all  intermediate 
points  r . 
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a 


2 


FIG.  3:  GFOMTRY  OF  II IF.  lTF.RATIYT  SOLUTION 

FOR  AN1  UNIFORM  ISOTOPIC  WALK 


Also  note,  that  since  the  one  step  probability  identified  above  is  properly  normalized, 


/ d£  I 6 [ 2 -ot2  ] = 1 


(C-4) 


as  one  can  easily  ascertain  by  explicit  integration,  the  normalization  condition 

j dr  G^r , (n- 1)  x | 0 , G ) - 1 (C-5) 

implies  via  (C-3), 

/ dp  G(p ,nx | 0 , 0)  - 1 
as  expected. 

The  very  reasonable  physical  interpretation  that  we  have  just  given  equation  (C  -3), 
suggests  the  possibility  of  guessing  at  the  results  to  be  expected  for  cases  that  are  some- 
what more  complex  than  the  simple  uniform  isotopic  case  analyzed  so  far.  Thus,  to  lift 
the  restriction  of  equation  (B-8),  we  allow  a,  (P)j  j to  be  any  arbitrary,  properly  normal- 
ized set  of  functions  and  write  for  the  probability  density  of  arrival  at  p in  n steps. 
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FIG.  4:  GHOMHTRY  OF  11 II • ITFRATTYT.  SOLUTION 

FOR  A NON- UNIFORM  ISOTOPIC  WALK 

That  (C-7)  is,  in  fact,  the  correct  answer  can  be  checked,  of  course,  with  the  tech- 
nique that  led  to  (C-3),  but  in  view  of  the  transparent  physical  content  of  (C-7)  that  is 
hardly  warranted 

Needless  to  say,  the  few  cases  mentioned  above  as  illustrative  instances  do  not  exhaust 
the  set  of  configurations  that  can  be  handled  by  the  iterative  technique,  but  we  shall  leave  it 
to  the  patience  of  the  interested  reader  the  job  of  looking  through  them.  The  one  configura- 
tion, however,  that  involves  relaxing-  assumption  (B-7)  as  well  as  the  ensuing  discussion 
are  considered  sufficiently  delicate  to  warrant  a separate  later  section.  As  to  the  manner 
in  which  actual  numbers  can  oe  obtained  from  the  rather  formal  results  of  this  section,  we 
shall  have  little  to  say  at  this  time,  but  note  that  quite  powerful  Monte  Carlo  methods  are 
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available  for  handling  exactly  this  kind  of  iterative  equations.  Correspondingly,  when 
needed,  solutions  for  any  value  of  n can  be  obtained  at  much  less  expense  than  experi- 
enced in  the  previous  section. 

D.  THE  CASE  OF  THE  CORRELATED  WALK 

For  all  their  value,  the  results  that  have  been  obtained  so  far  suffer  from  the  overall 
constraint  embodied  in  (B-7).  We  shall  now  show  that  the  main  features  of  what  under- 
standing has  been  gained  for  the  case  of  the  isotopic  walk  survive  the  introduction  of 
correlations. 

Consider  again, 

ai  * a ; vi  = v ; i = 1,2,  . . . ,n 
and  make  the  following  assumption  of  nonisotopy: 

a(Oj,02> • • • = 


to  replace  (B-7).  Here 


3-1 

represent  the  conditional  probability  distributions  for  sampling  0^  given  that  the  walking 
unit  is  at  the  point  in  the  ocean  of  coordinates, 

i - 1 1-1 

xi=  t a cos  G.  ; y.  = 2 a sin  0, 

3-1  J 1 j-1  J 

along  its  path.  Correspondingly, 


i-1 

I a cos  0 . 

3-1  J 


(D-l) 
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G ( p , n t | 0 , 0 j = fde,a}(Q  ) f 7 T d9  • 
J 1 J i=2 


a.  I & a 


i -J  i -1  \ 

X]  Ci  cos  0 . ; V a s in  e . * 

J = 1 J J = 1 J / 


(D-2) 


6[  23  a cos  © v ~ a ) «5(  I a sin  6.-0 

\k=l  k / \k= 1 k . 


or,  rather. 


+ r n /i'l  \ 

G(p,nt 1 0,0)  =J  7T  d0i  dy.  ax  (6i|xj,y.)  - £ a cos  6^  V 


(D-3) 


i-1 


n 


6^  - £ a sin  Q-j  6^  a cos  eR-o J ^ a sin  0,-p  J 

and  using  again  the  spectral  representation  of  the  Dirac  delta  function  together  with  the 
convolution  theorem  for  Fourier  transfer  ns,  we  obtain 


G(p,ni  | 0,0)  =J 


dr  2an(p-rir)  6[(p-r)-a  ] G^r, (n-l)x|0,0 ) 


where 


a 

n 


(p-r|  r) 


sin  0°  = l£ Hi 
n a 
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[t  then  follows  that  all  conceptual  as  well  as  calculational  advantages  that  go  along  with 
the  iterative  type  solutions  of  section  C,  remain  also  valid  for  the  correlated  configuration 
considered  here.  It  is  perhaps  worth  noticing  at  this  time  that  while  (D-l)  is  not  the  most 
general  correlation  along  the  path  that  can  be  envisioned,  it  should  account  for  the  majority 
of  realistic  situations  one  might  be  concerned  with.  Correspondingly,  we  shall  pursue  the 
question  of  correlations  no  further. 

E.  THE  CASE  OF  THE  BOUNDED  WALK 

In  many  realistic  cases  of  interest  to  us,  intelligence  information  is  available  con- 
cerning the  existence  of  a bounded  region  in  the  ocean  within  which  the  approaching  enemy 
vessel  is  most  likely  to  be  found.  Correspondingly,  the  study  of  the  arrival  probability 
density  in  the  presence  of  given  boundary  conditions  becomes  relevant.  Within  our  frame- 
work, the  question  of  boundary  conditions  can  be  treated  with  relative  ease  in  the  same 
spirit  the  question  of  correlations  along  the  path  of  the  walk  has  been,  i.e.,  by  translating 
the  physical  restrictions  imposed  therein  into  geometrical  constraints  for  (A-5). 

To  illustrate  the  basic  concepts  invo’ved,  consider  the  following  simple  boundary 
conditions, 

0 < x ,y  < L (H-l) 

corresponding  to  a square  shaped  limitation  upon  the  walk  of  the  target.  To  ensure  that 
(E-l)  is  respected,  it  shall  suffice  to  limit  the  set  of  paths  on  which  the  integration  in  (A-5) 
is  performed  to  include  only  those  that  independently  satisfy  (E-l).  Therefore,  if  the  case 
of  the  uniform  isotopic  walk  with  constant  speed  and  step  is  chosen  again, 


GBC(P,nx|b,0)  = 


'l\n  r n / 1 

— ) / 7T  dO  o(l  - la 

M 1 1 \ j-l 


cos  8 j 1 * 


n I a 


Of  I a cos  0.  1 0[  £ a sin  0.  I 01 i 

'yM  J / J/ 


I a sin  ()■ 
i=l  j 


(E-2) 


7n  \ Yn 

o(  T a cos  O.-u  A Z a sm  0.-p 
\j=l  J / \j=l  J 
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We  can  now  attempt  to  reduce  (E-2)  t'  an  iterative  solution,  employing  the  techniques 
developed  in  Section  C.  In  fact,  using  the  spectral  representation  for  the  step  function. 


e(x) 


-l 

2ri 


f 


- tJo 


doj 


ioix 
e 

ie 


one  can  show  that: 


GBc(p.nTl°.°)  * e(L-o)  0(0)  e(L-p)  e(p)  * 

J dr  ~ 6[(p-r)  -a  } G0C^r, (n-1) t|0,0 ^ 


(E-3) 


which  displays  Che  same  intuitively  simple  content  that  its  predecessors  have. 

The  attentive  reader  will  have  noticed  the  little  swindle  that  has  been  perpetrated  in 
the  definition  of  CBC ( P , n t | ft , Q ) * Specifically,  since, 

Kn  £ J dp  GBC(p,nT|ft,0)  * f dw  a(i:) 


where  the  right  hand  side  integral  extends  only  over  those  u paths  that  satisfy  (E-l),  the 
probability  density  is  not  normalized  to  unity.  In  fact,  one  can  easily  argue  that 


To  compensate  for  the  swindle,  let  us  define, 


vv 


' B C (■  p ’ n T i r0 


0) 


(E-i) 


ami  use  t!:e  symmetry  property 


G „r(p  » n'i  ! , 0 ) = GBC(r0,nx|  P,0) 


‘BC 


(K-5) 
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-\^-rrn 


to  write 


W 


J dp  Pp^ ( rQ > nx | p , 0 ) 


Hence,  in  virtue  of  (E-3) 

K0C?0)  = 1 


Kn(rc)  = 0(L-xo)  0(xo)  6(L-y0)  e(yQ)  * 


] d?  I 6[(r0-f)-o  ] Kn _x(r) 


an.  I for  each  value  of  n , the  properly  normalized  probability  density  is  obtained  dividing 
through  the  appropriate  solution  of  (E-3)  by  the  corresponding  solution  of  (E-7)  taken  at 


V°  ’ 


GRr(p,nx |ft,0)  = —— 


GRp(p,nx  | , 0) 


Kn(0) 


(E-3) 


' * 
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PARI  n.  DETECTION  PROBABILITIES 

Having’  modifie  1 the  Koopman  approach  to  account  for  non-uniform  probabilities  of 
arrival,  we  are  now  prepared  to  address  in  some  detail  the  problem  of  target  detection 
within  the  framework  of  a search  and  detection  game.  As  the  reader  will  recall,  the 
main  point  of  our  proposal  is  to  allow  the  target  to  arrive  in  the  vicinity  of  the  searcher 
along  any  one  of  a family  of  possible  paths  parameterized  by  the  multidimensional  set  uj  , 

F(r;  a;)  = 0 

Each  path  within  F is  subsequently  assigned  a given  probability  a (w)  of  having  been 
chosen  by  the  approaching  target,  and 

J dua(o)  = 1 
F 

Hence,  if  p(o),  t)  represents  the  probability  of  no  detection  through  time  T lor  a tar- 
get approaching  the  searcher  along  path  u>  , we  write  for  the  probability  of  no  det  ection 
through  time  t , 

p(t)  = f daia(w)p((u,  t) 

F 

A.  THE  DETECTION  ROBABILITY  FOR  A SPECIFIED  SEARCH  PATTERN 
Let  the  searcher  proceed  along  the  specified  search  pattern, 

rs  = rs(t)  . (AT) 


Then,  as  indicated  in  Koopman' s book,  the  probability  of  no  detection  through  time  t is 
obtained  by  solving  the  following  differential  equation: 


nfw.tl  v Tr . (m  tl 

i v v'  j — j i i “ 'j-  V — ’ j — j 


- r>i  i 

S v 1 


subject  to  tiie  initial  condition 

PO,0)  = p0  . 


Here,  and  everywhere  else,  (w  , t)  represents  the  position  of  the  target  along  its 
trajectory  at  time  t . 
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To  proceed,  we  take  Pq  = 1 , write  equation  (A -2)  as, 


P(w»t)  * 1 - J dr  Y [r^(oj , t)  - ?s(t)]  P(oj,t) 


(A -3) 


and  Newman  iterate  (A -3)  into  the  form, 

p((D,t)  = 1 + it,  (-l)n/  dx^  J 1 dx?  j 

n=l 


rT 

I n-1  , * 

dTn 


0 


0 


0 


(A -4) 


t[Y(?t(u>,t1)-?s(x1))  Y(?T(w,T2)-rg(T2))  . • -rsCxn))| 


In  equation  (A -4)  we  employ  the  symbolic  notation  T to  represent  the  time  ordering 
of  the  product  of  y ' s within  the  square  brackets,  i.e. , 

T «{>  (0  < t < T < . . . < Tn  < O 

We  make  use  of  the  symmetry  of  the  time-ordered  expression  upon  interchange  of 
any  of  the  coordinates  f 1 , x , • • • , Tn  . to  symmetrize  the  integration  interval  with 

respect  to  the  n indices.  For  n"2  we  recognize 


f dr1  fl  d^2  l|Yp  (a.,T1)-?sCT1))Y(Va,»T2)’rs(T2))J 


0 0 


t*  f 

f dT2  J 2 d‘l  T|Y(VW'Tl)'rs(Tl))Y(?T(a,’T2)'^Cl2))l 
0 0 


I /‘d’i  /*  *2  tHV-  ))t 


} 


J 


\ 

j 
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For  general  n we  can  similarly  carry  out  the  nl  permutations  cf  the  n indices  and 
extend  the  Integration  region  over  the  n-dimensional  cube  from  zero  to  t . Each  of  the 
n!  permuted  time -ordered  regions  contributes  equally,  and  we  therefore  write 


= i + / dix  / dT?  ...  r1  dtn 

n=l  n 0 Tl 


0 0 


t[y(?t(“»Ti)-?s(Ti))  • 

Next,  we  employ  the  identity 

YpT(w,T^rsCT.))  - /dr.  y(r.-rs(xij)  6 (r^.t^r.) 
to  write  for  the  probability  of  no  detection, 


,n 


OO  r -^1  **  ^ ^ 

pCoi.t)  = 1 + £ /*  drA  drx  | dx2d^2  ...  J1  dtd^  * 


n=l 


0 


" f 
0 


0 


(A -5) 


tKWt1>)  ^VVV)  6(V“-Tn)'?n)| 


Hence,  the  actual  probability  of  no  detection, 


p(t)  * f dot  a(oj)  p(u>,t) 


bicomes: 
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« ..  ,-.n  rt  t 

P(t)  = 1+1  J dTl  dxi  / dx2  dx2  ‘ ' 

n=l  ‘0  n 


/ 


dx^  dr 
n n 


T|Y(vrsCTl1)  Y(WT2>)  "•  Y(VYs(Tn>) 


(A -6) 


Jduta(tu)  «(r  (w.xp-rj)  fi(rT(u,T2)-r2)  . ..  6(rT^>Tn>rn)]  > 


where,  of  course,  in  obtaining  (A-6)  we  have  had  to  assume  that  the  infinite  sum  over 
n in  equation  (A -5)  lias  the  necessary  convergence  properties  to  justify  our  piecewise 

application  of  the  integral  operator  J dw  a (u>) 

In  equation  (A-6),  the  w integration  extends  only  over  the  family  of  solutions  to  the 
simultaneous  set  of  equations, 

rT(u>  ,x  ) - ri  (i  = 1,2,.  . . ,n) 


as  ensured  by  the  Dirac  delta  functions  appearing  there.  Within  this  family,  however, 
the  probability  corresponding  to  any  given  path  is  exprassable  as  the  product, 

n 

a (oj)  = 7 T (u>) 

i = l 


where  a.  represents  the  probability  of  selecting  from  all  trajectories  connecting 
u j -path  that  is  contained  therein 


(r^,  x^)  the  particular  one  that  coincides  with  the  segment  of  the 
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If  the  Markovian  type  assumption  is  made,  whereby  (w)  depend  only  on  the  set  uk 

parameterizing  the  family  of  curves  connecting  (r^_^  , t.  ^)to  (r^>  T^)  * ^ can  ^ 
easily  shown  that, 

/ du  a(w)  ... 


~ 7 7 g(t - , t.  r.  , t.  ) 

' 1 i i-l  i-l  / 
x=l 

*>  / 
for  Tq  the  initial  target  position,  Tq  = 0 , and  where  Gyr^,  r^  ^ j 

are  the  probability  densities  of  arrival  discussed  in  Part  I of  this  paper. 

Consequently, 

°°  11 

p(t)=l+  I # dTjdr.  v(?r..?s(ti))  ‘ 


0 i-l 

It  is  perhaps  worth  noticing  that  the  disappearance  of  the  symbolic  time  ordered 
product  in  equation  (A -7)  is  only  formal  in  that  each  ^(r|>Ti}r|  i»T- 

implicitly  carries  with  it  a causal  step  function,  0 / r . - 1 . 1 \ ensuring  that  all  propa- 
gation takes  place  forward  in  time . ' x 1 ' 

One  way  to  render  equation  (A -7)  useful,  would  be  to  define: 

(-+  ■+  \ £ -*•  -y 

r ,t  r ,r  ) = / dx  dr.  — / dx^dr  * 

n+1  n+1  0 0/  J 1 1 J n n 
0 0 

(A -8) 


G (r  1 , t 
V n+1  n+1 


i'To)  it  *(V^Mc(ViFi-i-Ti.i)  ’ 
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and  correspondingly  write  for  the  probability  of  no  detection, 


uu  H "t 

pm  • 1 - r 7 tyty  ^Yh'V)  Vl(  vTnlvTo) 

n=T  0 


(A -9) 


Although  not  necessarily  transparent,  there  is  virtue  in  having  written  p(t)  in  the 
form  of  equation  (A-9).  The  advantage  is  contained  in  the  observation  that  thus  defined 
satisfy  an  iterative  type  equation  of  the  form 


♦n(Vl'VllVo)  ' / dTnrfn  '(Vl’VjV'n) 


(A -10) 


with, 


♦opi-hl'o-'o)  VTu)  • 


(A -10') 


In  fact,  if  we  choose  to  identify  ^n(^n+l  ’ T]t-t  1 | ro  ' Tq)  witl1  tlie  P^k^bility  tar" 

get  arrive1  at  ^n  + l,Tn  + l)  the  presence  of  an  n~uple  attempt  at  detection,  equation 

(A-10)  obtains  the  following  simple  physical  interpretation  (see  figure  5):  the  probability 
of  arrival  at  f rn+  ^ , t^  + ^ ) in  the  presence  of  an  n-uple  attempt  at  detection, 

(-+  | -+  \ -y 

rn  + l,Tn+l  I r o » T o / ,-)  i>  *s  0l*ained  by  multiplying  the  probability  of  arrival  at 

some  point  / i T ) in  the  presence  of  an  (n-l)-uple  attempt  at  detection, 

'f'  i / 1’  , t rn  iin|  dr  » with  the  instantaneous  probability  of  detection  Y ( r - r (t  ))  dT 
nxyrifiUO/n  ^vn^n/n 

by  one  glimpse  at  , propagating  the  product  thus  described  from  ^ rR  » T n ) to 

(*n+l,Tn+l)  viathe  "propagator"  G(?n  + 1 » Tn  + 1 | ’ Tn  )d^n+l  and  then  Integrating 

over  all  the  allowed  Intermediary  points  (^n’  Tn  J 


3 
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FIG.  5:  GEOMETRY  OF  THE  ITERATIVE  SOLUTION  FOR  DETECTION 
ALONG  A SPECIFIED  SEARCH  PATTERN 


It  is  quite  important  to  recognize  the  crucial  role  that  the  propagator  plays  in  the 
formalism  discussed  above;  it  is  precisely  the  presence  of  G^rn+1,Tn+^|rn,Tnj 

in  equation  (A -10)  that  properly  accounts  for  whatever  correlations  there  exist  along  the 
path  of  the  searcher,  in  that  the  propagator  is  responsible  for  the  target  arriving  at  the 
sight  of  the  (n+l)th  attempt  at  detection  with  full  memory  of  the  point  ^ rn , tn  j where 

the  nth  attempt  at  detection  took  place.  In  fact,  if  one  were  to  arbitrarily  break  the  corre- 
lation by  replacing  G(rn+1,Tn  + 1 rn,tn)  in  equation  (A -10)  with  G (^n+1 » Tn  + 1 | » To  ) 

the  entire  framework  would  automatically  collapse  into  a Kooprnan-type  result  for  p(t)  . 
Indeed,  upon  defining 


=/  dT„d?„  Vi(vTn|vTo ) 


(A -11) 


and  performing  the  indicated  replacement,  equation  (A-10)  becomes 


0n(t)  *0^(10  J drdr  y(t-?s(t))  C'(r,x|  ?0,t0  j 


(A -12) 
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and  correspondingly, 


P(tj  =1+1  ^ ,Ct)  • 

n=l  n!  n'l 

But,  in  virtue  of  equation  (A -12) 


(A -13) 


<pn(t) 


J dt  dr  y^r-rs(x)^  G (r , t J r0 , tq) 


n + 1 


(A -14) 


so  that 


’ 1 * „fi  ^ 


J drdr  Y(r-rg(T)J  g(^t  |r0,T0) 


(A -15) 


or,  rather, 


p(t)  = e 


f ^ ■+  / ^ , 

/ dtdr  y(r-r  (t) 
J0  ' 

l)  c(r,T|r0,Tf) 

(A  - 1 6) 


to  be  compared  with  equation  (1-3).  If  one  is  further  willing  to  assume  that  G^r , t ' to) 

represents  uniform  arrival  within  an  ocean  region  of  area  A,  equation  (A-16)  simply 
reduces  to  (equation  (1-4))  as  expected. 

Returning  to  the  general  formalism  as  comprised  of  equations  (A -9)  and  (A -10),  we 
finally  note  that  it  displays  the  appropriate  str  ucture  to  render  the  Monte  Carlo  numerical 
methods  mentioned  in  Part  I equally  useful  to  the  detect' on  problem  studied  here. 

B.  GLOBAL  DETECTION  PROBABILITIES 

In  the  foregoing  discussion,  the  path  of  the  searcher  has  been  arbitrarily  specified 
and  the  principle  of  biasing  the  set  of  possible  paths  in  accordance  with  a(cu)  used  only  in 
the  description  of  target  motion.  The  distinct  advantage  of  doing  so  lies  in  the  possibility 
of  using  the  formalism  thus  developed  for  analyzing  search  pattern  optimization  problems. 
Sometimes,  however,  there  is  value  in  calculating  global  probabilities  of  detection,  where- 
by both  the  searcher  and  the  target  motion  are  treated  via  the  biasing  principl  . When  such 
is  the  case,  we  shall  use 

dpCo,w;t)  = - p(a,w;t)  y[rT(u,t)  - rs(o,t)  ] dt 
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(B-i) 


6 


to  replace  equation  (A -2),  where  o and  u>  represent  the  family  of  trajectory  parameters 
for  the  searcher  and  target,  respectively.  Therefore, 


p ( o , u , t ) 


0 


-?s 


and  Newman  iterating  (B-2), 


p(o,u;t)  = 1 + T 


Izll 


1 ~T  f ^ i f dT 

n=i  n-  7 1 J n 


(B-2) 


(B-3) 


Just  as  before,  we  now  use 

Y(rT(oJ,T.)-?s(o,T.)  j - 

fd*idpi  £(TT(ui>Ti)~*i)  6(rs(a’Ji^^i) 

to  write  equation  (B-3)  in  the  form: 

oo  n 

_ i-i)  ,t  __  t 

p(CT,o>;t)  = 1 + In!  j (dT1dr1dp1)  ...  f (dx  drndpn)  * 

n“1  -Q  JQ  n 


t[y(Vpi)6(VUiTiv?i)  6{Va’Ti)-pi)  ••• 


/n_  A\ 
\11 


Y(rn- 


n-pn)  d(rT(a,,Tn)-rn)6(?s(a,Tn)-Pn)]  . 
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Hence, 


p^)  = 1 + ’ J (d-r1dr1dp1)  ...  (dTndrndpn)  * 

0 

Tlrfrj-Pj)  ...  Y(rn-pn)  J dusdcs  A(o,u)  * 
6(rT(«,T1>5^<(rs(o.T1)-pjj  ...  i(Vw>Tn)‘V)fs(0-Tn: 

where  A(a,aj)  stands  for  the  joint  probability  that  the  searcher  follow  path  o 
target  path  u . 

We  shall  assume  for  simplicity  that: 

ACa,u>)  = AjCo)  A (w) 

2 

which  is  tantamount  to  ignoring  the  dynamics  of  the  search  and  detection  game 
Correspondingly,  s 


SzXH  * -*  -*■  - 


pw  - 1 + ; ~n:  / cdT^dPi) ...  / (dxnd?n 


d»n>  * 


A ^Vdi>  Gi^i’Ti!di-l’Ti-l)  Gs(pi,Ti'i’i-l'Ti-i  ) 


(B~  5) 


ana  the 


03-6) 


(B-7) 
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Again,  we  i itroduce: 


*n(VrW'Tntl'ro’pO’’o)  : 


f tdTldi?ldPl)  •••  (dTndW  SlVl'Vl'V'J'ill'Vl'V'.) 


7T  YC?.-pJ  G (r  ,x  |r.  t.  ) G / p. , . |p.  ,t.  ) 

i=1  i 1 A i l1  l-l  1-1/  s\wi  l i-i  i-i  / 


to  write. 


oo  t 

— ^ f — ► — V-  4-  ->  ~>”f-  ->--V 

p (t)  , 1 + L_J_  j (dTndrndpn)  Y(rn-pn)  V1(rn’Pn'Tnlr0*Po'To)  (B'9) 

and  notice  that  satisfy  the  iterative  equation,  (see  figure  6), 
n 

tn(Vl’Vl,Tn+lirO’P0’To)1 

r -+  ->•  /-*•  \ ,-*■  ■+•  (B-10) 

/ (dr  dr  dp  ) G ( P ,t  i p ,t  G (r  ,x  |r  ,t  ) * 

J n n MnJ  s\  n+1  n+1 1 n n'  T\  n+1  n+1  n n ’ 


y(r  -p  ) (r  ,p  ,t  |r  p x ) 
\ n n / n-l\n  n n'O  0 0/ 


1 


r 


i 

i 


FIG-  6:  GEOMETRY  OF  THE  ITERATIVE  SOLUTION 
FOR  GLOBAL  DETECTION  PROBABILITIES 


Naturally,  equation  (B-10)  contains  physical  information  similar  to  that  identified  in 
equation  (A -TO)  and  therefore  most  of  the  discussion  following  equation  (A -10)  should  apply 
here  too.  Also,  by  the  nature  of  (B-10),  some  generalization  of  the  Monte  Carlo  methods 
suggested  there  for  the  solution  to  the  specified  search  pattern  problem  ought  to  be  useful 
in  finding  global  probabilities  of  detection. 
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DISCUSSION 


Motivated  by  the  need  for  developing  a framework  of  larger  breadth  and  flexibility  than 
the  Koopman  one,  we  have  succeeded  in  delineating  a train  of  thought  that  holds  the  promise 
of  successfully  accounting  for  the  modern  sophistication  of  the  search  and  detection  game. 

Thus,  a principle  of  biasing  has  been  employed,  whereby  given  probabilities  are 
assigned  to  all  possible  paths  of  a military  unit,  and  the  probability  of  arrival  of  that  unit 
at  any  point  in  the  ocean  expressed  in  terms  of  them.  It  then  becomes  possible,  as  shown 
in  Section  C,  to  formulate  the  result  in  terms  of  iterative  solutions  whose  physical  content 
is  interestingly  simple  and  whose  calculational  advantages  are  many  fold.  In  fact,  a full 
solution  can  be  obtained  numerically,  if  the  input  probabilities  per  path  are  given.  Whether 
they  reflect  our  total  ignorance  concerning  the  unit's  motion  or  the  h*ghly  correlated  nature 
of  the  path,  a choice  of  such  inputs  is  always  possible,  in  that  the  inputs  merely  represent 
our  actual  data  base  ; nd  not  some  restrictive  model  assumption  that  might  not  correspond 
to  the  reality  we  want  to  describe.  The  flexibility  that  this  imparts  to  the  formalism  is  such 
that  one  might  approach  within  the  same  framework  problems  that  have  varying  amounts  of 
input  information. 

Let  us  also  note,  that  by  the  very  nature  of  the  iterative  solutions  we  propose,  a con- 
nection can  be  established  between  the  search  and  detection  problem  on  one  hand  and  the 
well  studied  mathematical  theory  of  Green's  function  on  the  other,  thus  opening  a door  for 
the  influx  ot  a whole  family  of  numerical  methodologies  that  have  proven  useful  elsewhere. 
Although  much  might  still  be  said  about  the  hidden  or  obvious  dimensions  of  this  frame- 
work, we  stop  our  discussion  here  for  we  wish  merely  to  enlighten  our  reader's  judgement, 
not  to  influence  it. 


